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Abstract. 

Near-infrared surveys of high-mass star-forming regions start to shed light 
onto their stellar content. A particular class of objects found in these regions, 
the so-called massive Young Stellar Objects (YSOs) are surrounded by dense 
circumstellar material. Several near- and mid-infrared diagnostic tools are used 
to infer the physical characteristics and geometry of this circumstellar matter. 
Near-infrared hydrogen emission lines provide evidence for a disk-wind. The 
profiles of the first overtone of the CO band-heads, originating in the inner 10 
AU from the central star, are well fitted assuming a keplerian rotating disk. 
The mid-infrared spectral energy distribution requires the presence of a more 
extended envelope containing dust at a temperature of about 200 K. CRIRES 
observations of CO fundamental absorption lines confirm the presence of a cold 
envelope. We discuss the evolutionary status of these objects. 



Introduction 



The observational study of high-mass star-forming regions is hampered by the 
fact that they are so unique, and thus, on average distant. Additionally, the 
actual formation happens deeply obscured inside molecular clouds behind hun- 
dreds of magnitudes of visual extinction. The actual formation sites are only 
observable at radio and (sub)mm wavelengths. These observations, however, 
only provide indirect information about the forming stars. The Lyman contin- 
uum photon rate derived from the radio flux or the integrated infrared luminosity 
provides an estimate on the properties of the massive star(s). Only after the 
stars are formed and start to clear out their environment, they become observ- 
able at shorter wavelengths. The near- infrared window is ideally suited to study 
the massive stars in this evolutionary phase. The photosphere and close circum- 
stellar environment of the young massive stars can be detected as the extinction 
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is modest and the thermal emission of the surrounding molecular cloud is not 
dominant. 

Recent imaging and spectroscopic surveys have revealed the stellar con- 



tent of high-mass star-forming regions (lHanson et al.ll2002l: iKendall et al.ll2003l: 



Alvarez et ap2004l : lBik]l2004 : Blum et al.ll2004l : lBik et al]l2005l . l2006l : lKaper et al" 



20071 ). Several types of objects were classified in these complex reg ions. A large 



number of OB stars have been identified by their photospheric lines ([Watson Hanson 
19971 : lHanson et al.ll200l : iBik et aLll2005r ) and their properties do 



not seem to 

differ from OB main sequence field stars observed in the optical. Another class of 
high-mass objects identified are the massive Young Stellar Objects. Their near- 
infrared colors indicate an infrared excess, and their near- infrared spectra are 
dominated by em ission lines from ionized (hydrogen, h e lium) as well as neutral 



species, e.g. CO (|Kendall et al.ll2003l : iBlum et alJ[20M : IBik et al.ll2006l l 



These spectroscopic features as well as the infrared excess can be used to 
derive the physical parameters of the cirumstellar material and try to understand 
the evolutionary status of these objects. The detection and study of disks around 
young massive stars that are remnants of the f ormation process is key t o our 
understanding of the formation of massive stars ( Zinnecker &: Yorke 20071 ). 



2. Near- and Mid-Infrared diagnostic tools 

In this section the different diagnostic tools available in the near and mid- infrared 
window will be discussed. Line ratios of ionized hydrogen lines are used to trace 
the physical conditions and geometry of the ionized circumstellar material. The 
line profiles of the CO first overtone emission provide kinematic information 
about the dense, neutral gas in the inner 5-10 AU. Extended envelopes, located 
further away from the central star can be detected at longer wavelengths. 

2.1. Ionized gas 

One of the spectral characteristics of the massive YSOs is that they show lines 
emitted by the recombination of hydrogen and helium. In the K-band, not 
only Br7, but also weaker lines of the hydroge n Pfund series ( Pf23 - Pf 33) are 
observed, as well as Hel lines in some cases (Bik et al.l l2006). The Hel lines 



indicate the presence of a hot io nizing source with a spectral type of at least 
late - mid O ([Hanson et al.ll2002l ). 



The hydrogen lines can be used to obtain an estimate on the density of the 
emitting gas. The Pfund lines are high-density tracers which are only emitted in 
very dense gas (N e up to 10 9 cm" 3 ). These lines likely originate in a circumstellar 
disk. This is, moreover, supported by the very broad line profiles (300 - 500 
km s _1 ) or by double peaked profiles in a few objects. Br7, on the other hand, 
can be emitted in a large range of environments with different densities. It can 
be associated to diffuse HII regions, stellar winds, disk-winds, or very dense 
gaseous disks around Be stars. Also here, the observed spectral profile can help 
to discriminate between possible scenarios. The observed velocities of Br7 are 
in the range of 100 - 200 km s . This excludes the possibility that the lines are 
of nebular origin. The recombination lines in an HII regions have a broadening 
similar to the local sound speed which is on the order of 20 km s" 1 for ionized 
gas. Also, the origin in a stellar wind can be excluded as the velocity profiles 
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Figure 1 . left 2 pannels: L-band spectra taken with ISAAC at the VLT of 
12 massive YSOs. The broad absorption are molecular solid-state transition 
features caused by material in the diffuse ISM or molecular clouds. These 
absorption features can be used to derive an accurate line-of-sight extinc- 
tion value. The emission lines are mainly hydrogen emission lines. The lines 
needed to construct the diagnostic diagramme are Bra(4.05 /an), Pf7(3.754 
fim) and Hul4 (4.02 (im). right: The L-band diagramme as described in 
the text. The gray circles repres ent the loci of t h e mas sive stars with known 
cirumstellar geometry for which iLenorzer et al.l (|2002f ) constructed the dia- 
gramme: LBVs, B[e] and Be stars. Overplotted are the massive YSOs. 



of stellar winds show very broad lines (500 - 2000 km s 1 ). Double peaked Br7 



emiss ion in some massive YSOs clearly favors the disk-scenario (|Blum et al 



2004). However, not all the Br7 lines show this profile and a contribution of the 
disk-wind cannot be excluded for these. The observed velocity of the Br7 line 
is consistent with the disk-wind scenario. 

A disk-wind is caused by the interaction between the UV photons and the 
ionized upper layer of the disk. This result s in a radiation d riven, outflowing 
disk-wind with velocities up to 200 km s~ 1 ( Drew et al. 1998). In a disk-wind, 
the average density of the emitting ionized gas is lower than in a geometri- 
cally thin circumstellar disk. This will affect the relative strength of the hy- 
drogen lines. While in low density regions (e.g. in an HII region), the relative 
strength of hydrogen lin es follows the optically thin (case B) approximation 
( Storey Sz Hummer 1995j); in high density regions, where the gas is optically 
thick, the relative strength of the lines is determined by the emitting surface of 
the gas. 



Lenorzer et al. (2002) constructed a diagnostic diagram using 3 specific hy- 



drogen lines in the L-band (Bra, Hul4 and Pf7) to demonstrate this method. 
The ratio of Hul4 over Bra spans a range of 2 orders of magnitude (from ~ 0.02 
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Figur e 2. Observed CO bandhead profiles from 2 massive YSOs (|Bik fc Thil 
2004). The top spectra are the observed spectra, the bottom spectra are 
the best fitting models. These models indicate that the CO is located in a 
circumstellar disk within 5 - f AU from the central star. 



in case of an optically thin HII region to unity in the case of an optically thick 
emitting medium) . In order to calibrate the diagr am and to relate it to different 
circumstellar environments, Lenorzer et al. (|2002l ) used the observed line ratios 



of massive stars with a known circumstellar geometry (Fig lb). Three different 
types of objects were chosen; Luminous Blue Variables (LBVs), evolved massive 
stars with a very strong stellar wind, B[e] stars with an equatorially flattened 
dense stellar wind or disk and the Be stars with an optically thick, geometrically 
thin ionized disk. The location of these respective objects is plotted in Fig lb as 
circles. The Be stars fall as expected close to the optically thick location while 
the LBV stars are located more in the direction of the optically thin material. 

We observed a sample of 12 massive YSOs in the L-band with ISAAC at the 
VLT (Fig. la) and measure the line ratios to compare them with the location 
of the well studied massive stars in the diagram. In Fig lb the massive YSOs 
which have all the 3 hydrogen lines in their spectrum are overplotted on the 
diagram. 

The location of the massive YSOs overlaps that of the B[e] and Be stars. 
First this confirms that the ionized emitting region of the massive YSOs is, at 
least partially, confined into a disk. The span in location implies large differences 
from object to object for the averaged densities of the ionized circumstellar gas. 

A possible explanation for this, is that some objects are surrounded by a 
geometrically thin and confined disk (similar to Be disk) while others have disks 
that are more puffed up as expected when the disk is under the process of being 
photo-evaporated. 



2.2. CO first overtone emission 

Additional to the ionized lines, lines emitted by neutral species are observed in 
the spectra of massive YSOs. Among the most prominent features in the near- 
infrared are the CO first-overtone emissi on lines. These bands are observed in 
about 25% of all the massive YSOs (e.g IChandler et al.l fl993l . Il995l : iBik & Thi 
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Figure 3. Spectral energy distributions of the candidate massive YSOs. The 
data points are fit by a simple model. Panels a,b,c and e show a model fit 
with a 2 temperature black body, the data points of panel d are fitted with a 
powerlaw and the fit in panel f is a single temperature black body. 



20041 : iBlum et al.ll2004l L The CO bandheads are emitted in neutral material 
with temperatures between 2000 and 5000 K and densities of about 10 10 cm -3 . 
This hot neutral gas is located in the inner, dust-free regions of the circumstellar 
environment, relatively close to the star. 

The kinematic information provided by high resolution spectra of the band- 
heads can be used to constrain the location of the emitting material in the cir- 
cumstellar matter. Moreover, intrinsically, the blue side of the bandhead is very 
steep and different velocity profiles give rise to a different shape of the bandhead. 
Keplerian disks create a blue wing on top of the profile ( Fig 2). Thi s blue wing 
cannot be explain ed by a gaussian velocity distribution. iBik Sz Thi ( 20041 ) and 
IBlum et al. ( 20041 ) succesfully fit the high-resolution CO spectra with a emission 
profile from a rotating disk, assuming optically thin emission. 

This model also shows that the CO is emitted very close to the central star, 
in the inner 5-10 AU of the circumstellar environment. At these distances from 
the central star, the CO molecules need to be shielded from the UV photons in 
order not to be photo-dissociated. Therefore, this gas cannot be located at the 
surface of the disk, that would be ionized by the stellar radiation. The optical 
thickness rises steeply in the radial direction, shielding the midplane of the disk 
from the UV radiation, allowing the material to become neutral and the CO 
molecules to survive. 



2.3. Spectral energy distributions 

One of the main characteristics of massive YSOs is that they have an infrared 
excess. Their Spectral Energy Distribution (SED) is not dominated by pho- 
tospheric emission, but by the emission of circumstellar matter, mainly dust 
emission. The near-infrared emission lines, as discussed in the previous subsec- 
tions trace the ionized or hot neutral gas in the circumstellar environment. The 
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Figure 4. Absorption profiles of 3 12 CO lines observed with CRIRES at high 
spectral resolution (6 km s _1 ). Two absorption components are detected, the 
strongest component at +5.5 km s _1 is saturated, the weaker component at 
-5.5 km s" 1 is only detected in the R(4) to the R(0) lines, suggesting a very 
low temperature. 



infrared SED traces the warm material (down to a few 100 K) in the cirmstellar 
environment located further away from the central star. 

We have performed photometry of a sample of massive YSO covering the 
wavelength range from 1 to 20 fim using the ESO instruments SOFI, TIMM2 
and VISIR (see Fig. 3 for some exam ples). The observ ed SEDs are dereddened 
using extinction values determined bv lBik et al.l (|2006l ). The SEDs show a large 
range in spectral slopes, suggesting large differences in the properties of the 
circumstellar environment. They vary from blue slopes (Fig 3, panel d) to very 
red slopes (panels b,c and e). 

The SED of the bluest object in our sample (11097nr693) cannot be ex- 
plained by dust emission. The blue slope can only be produced by optically 
thick free-free emission from e.g. a gaseous disk like in Be stars. The SEP ca n 
be fitted with a powerlaw, with a slope typical for Be stars ( Waters et al.|[l99ll ). 
Surprisingly, the SED around 10 fim suggests the presence of silicate emission, 
not expected in a gaseous disk. This could be caused by a contribution of the 
high background from the surrounding HII region. 

The other objects have a much redder SED, and their SED can be explained 
by dust emission. Overplotted to the observed data in Fig. 3 is the best fit using 
a one or two temperature blackbody. The objects fitted with a 2 temperature 
blackbody show that a colder component (T ~ 200 K) tracing dust further 
away is needed to explain the shape of the SED. These objects also are spatially 
resolved at mid-infrared wavelengths and their observed size becomes larger with 
wavelength. For one of the reddest objects, 16164nr3636, the size of the 10 and 
20 /im emission extends up to 7500 AU away from the central star. Far-infrared 
observations of this source reveal a strong point source with a derived mass of 
about 2OOM0 (jKarnik et al.ll200lh . 



2.4. Absorption spectroscopy with CRIRES 

Another way to study the outer regions of the circumstellar material is the 
analysis of absorption lines of e.g. CO towards a bright backgro und. As part 
of the Science Verification observations ( Siebenmorgen et al. 2007) for CRIRES, 
one of the latest VLT instruments, we observed one massive YSO (IRAS 16164- 
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Figure 5. A schematic view of the different line- forming regions in the cir- 
cumstellar disk surrounding a massive YSO. The CO is located in the dense 
mid plane of the disk, while the hydrogen lines are emitted by the ionized 
surface layer and/or disk wind. 



5046) at 4.5 fim to detect the 12 CO lines of the fundamental transitions. The 
spectral resolution of these observations is R = 50,000 (6 km s _1 ). The 12 CO 
lines we observe are in absorption and show multiple components. In Fig. 4, 
the R(7), R(4) and R(0) lines of 12 CO lines are plotted. In all the spectra, 
a saturated absorption profile is present at -5.5 km s -1 . In the spectrum of 
the R(4) line a second absorption is present at +5.5 km s _1 . This line is not 
present in the R(7) spectrum but it appears even stronger in the R(0) line. This 
suggests that this line is absorbed in very cold mat erial of a few 1 K. The cold 
envelope detected in the far-infrared observations of lKarnik et al.l ( 200ll ) is likely 
the material that causes the saturated absorption profile of the 12 CO lines. The 
other component could be related to the source, but also to a cold molecular 
cloud in the line of sight. 



3. Summary and Evolutionary status 

Several different diagnostic tools are discussed to obtain information about the 
physical and kinematical properties of the circumstellar material. The different 
diagnostics trace different physical conditions and therefore different locations 
in the circumstellar environment of the massive YSOs. Based on the collection 
of observations a sketch can be drawn on how the circumstellar environment of 
the massive YSOs could look like. In Fig. 5 the different environments and their 
distance to the central star are sketched. Note that not all the objects have all 
the characteristics discussed here, some only show a subset of them. 

Stellar photons ionise the surface layer of the inner disk, which recombina- 
tion is observed in the hydrogen and Hel lines; the radiation pressure exerced 
onto the surface layer also drives matter away in the form of a disk wind. In 
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some objects the inner parts of the disk are so dense that the mid-plane can 
become neutral as it is shielded for the UV photons. In these regions (the inner 
5-10 AU) originate the CO bandhead emission showing evidence for keplerian 
rotation. Further out, the midplane of the disk cools even more and dust is 
able to survive giving rise to a red SED observed in the mid-infrared. A few 
objects show evidence for a cold dust component suggesting the presence of a 
more extended envelope. Most objects, however, are not detected in the mm 
regime, suggesting that they lack a large reservoir of very cold dust or gas. 

This illustrates the difference with the objects discussed by Zhang (2007) 
where large disks are detected in the (sub) mm regime. Still in the process of 
formation, these objects are too embedded to be detected in the near-infrared. 
Likely these large and cold outer regions of the circumstellar material are (being) 
photo-evaporated away. Photo-evaporation first leads to the disruption of the 
outer disk, while the inner disk can survive much longer (Hollenbach et al. 
1994). The inner disk is what is detected at near- and mid-infrared wavelengths. 
Objects for which the observations indicate the presence of an envelope might 
be the younger objects of our sample. Another important factor in the photo- 
evaporation process is the amount of UV photons emitted by the central star. 
This can drastically influence the dispersion time scales. A larger sample would 
be needed to disentangle these effects. 
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